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^ ■ Abstract 

o 



The recent extension of the standard model to include massive neutrinos in 



o , 

, the framework of noncommutative geometry and the spectral action principle 

' involves new scalar fields and their interactions with the usual complex scalar 

doublet. After ensuring that they bring no unphysical consequences, we address 
^ I the question of how these fields affect the physics predicted in Weinberg-Salam 

^ theory, particularly in the context of the Electroweak phase transition. Apply- 

ing the Dolan-Jackiw procedure, we calculate the finite temperature corrections, 
and find that the phase transition is first order. The new scalar interactions sig- 
nificantly improve the stability of the Electroweak Z string, through the "bag" 
phenomenon described by Watkins and Vachaspati. (Recently cosmic strings 
have climbed back into interest due to new evidence). Sourced by static em- 
bedded strings, an internal space analogy of Cartan's torsion is drawn, and a 
possible Higgs-force-like 'gravitational' effect of this non-propagating torsion on 
the fermion masses is described. We also check that the field generating the 
Majorana mass for the vr is non-zero in the physical vacuum. 



*Email: rmartins@math.ist.utl.pt, research supported by Fundagao para a Ciencias e a Tecnologia 
(FCT) including program POCI 2010/FEDER. 



1 



Contents 



1 Introduction 

1.1 Background information on noncommutative geometry 



2 Field content 

2.1 Introduction [sl 

2.2 Lorentzian [o] 

2.3 Euclidean EH 



3 Quantum and finite temperature corrections in </> theory [14 

3.1 Phase transitions [15 

3.2 The Coleman- Weinberg Mechanism Il6 

3.3 Finite temperature corrections 



4 Analysis of the effective potential: determining the order of the phase 



transition in two theories 23 



4.1 Analysis of the effective potentials [24 



4.1.1 Finite temperature corrections |27 



4.1.2 Critical temperature |2 



4.1.3 The Ehrenfest classification 29 



4.1.4 Type of extremum in the symmetric phase l30 



4.1.5 Summary |31 



5 Embedded defects and torsion 32 



5.1 Introduction to topological and embedded defects |32 



5.1.1 What are topological and embedded defects? |32 



5.1.2 Why do they form? |33 



5.1.3 How is torsion sourced by line defects? |34 



5.2 Electroweak Z string solution |34 



5.3 Stability of Electroweak strings |35 



5.4 The "Higgs force" and torsion gravity |37 



5.5 Vacuum manifold metric and string energetics |3; 



1 Introduction 



The noncommutative geometric formulation of the standard model coupled to gravity 
by Connes and Chamseddine (the spectral action principle |SAPj ) has been extended 
to include massive neutrinos jSMVp . [LNG], jSNM] [GM]. 

By definition, the standard model describes 3 massless chirality left neutrinos vl and 
precludes chirality right neutrinos vr. Further to the results of laboratory experiments 
it is nowadays consensual that neutrinos are massive particles and due to lack of 
evidence, the question of the existence of is open. Theoretical studies of particle 
physics now include appropriate generalisations of the standard model, but as the 
theory of Weinberg, Glashow and Salam has been shown to be very strong in its 
predictive power, the modifications are made as minimally as possible. 

In particular, the generalisation referred to above is done by the omission of the unax- 
iomatic assumption of S'°-reality and the natural inclusion of vr. However, it cannot 
be said that nothing else in the model may alter as a consequence; the most notable 
new feature being the appearance of additional scalar fields. These are note-worthy be- 
cause such a change can radically alter a theory: for example, it is normally considered 
undesirable to involve new fields, as they are likely to contradict experiment. 

Indeed in the first studies of such new fields |NGLj called leptoquarks, they cause colour 
symmetry breaking, which is of course absolutely unacceptable. More positively, in 
the extension being discussed here, some of the new fields are colour-charge-free, and 
might be as difficult to detect as vr, while the leptoquarks are either non-dynamical, or 
ruled out altogether by an additional constraint applied to the Dirac operator [GSN] . 

To add scalar fields to the theory is by no means bound to be a trivial alteration 
and given that the extension does not contradict experiment, let us instead ask the 
affirmative question of what predictions it may lead to. 

Recently there has been a lot of fresh interest in cosmic strings since two things have 
changed over the last few years. Cosmic strings are expected because of their oc- 
currence in many supersymmetry grand unified theories, but due to inconsistencies 
between predictions and the cosmic microwave background temperature anisotropics, 
it was necessary to explain why they and other topological defects should exist but 
do not. The favourite explanation is that they were "infiated away" beyond observ- 
able horizon. Topological defects were once thought to provide the seeds of structure 
formation, but these data ruled them out. The two new things that have happened 
|RCSj are (a) with recent advancements in theoretical physics, cosmic superstrings are 
expected to play the role of cosmic strings, and (b) more importantly, the data that 
ruled them out is being revised. Recent evidence implies that the parameters may be 
constrained sufficiently to allow agreement between predictions and measurements. 

At the level of Weinberg- Salam theory, there are no topologically stable defects as 
the vacuum manifold, the three-sphere, is simply connected. Nevertheless solitonic 

^The validity of the Majorana type terms in found in this article should be considered modulo the 
fact that the calculations were done in the Euclidean signature. 
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solutions there are, which are referred to as embedded strings and sphalerons. These 
embedded defects ahhough classically stable, have been shown to be unstable to small 
perturbations within a realistic range of parameters. Watkins and Vachaspati [BESj 
realised that an extension to the standard model by additional scalar fields forming 
bound states with the Electroweak Z string can improve its stability. It is clearly 
appropriate to ask whether the new fields in the extension we are studying can con- 
tribute to the stability in a similar vein. What we find is that the equations of motion 
precisely describe the scenario predicted in |BES] . 

The following questions are motivated: 

• What is the precise field content of the extended theory, and does this generali- 
sation have unphysical consequences? 

• Is the mass generating field for the v non-zero in the physical vacuum? 

• What effect if any do the new scalar fields have on the order of the Electroweak 
phase transition? In other words, do they ensure a gradual or a discrete change 
in the phase? 

• Do they have any stabilising effect on the embedded strings? 

• What is the role played by these fields in the context of internal space physics 
as described by Connes in |GM] ? That is, in the equivalence principle where the 
Higgs takes on the role of the gravitational field strength as a connection on the 
dimensional internal space. 

In the first section we establish the total scalar field content of the theory, explaining 
why there is no conflict with experiment, and show that the new scalar that generates 
the Majorana mass for is non-zero in the physical vacuum. We also discuss a Wick 
a rotation from the Euclidean signature to the Lorentzian for the case of the model 
with i^'O-dimension 0. Secondly, we introduce the Coleman- Weinberg |CW] mecha- 
nism for calculating the quantum corrections at zero temperature in </)^-theory, and 
then work through the details in the Dolan-Jackiw calculation of the finite tempera- 
ture corrections [DJ] in the same theory. The latter procedure is applied to the said 
noncommutative geometry extension of the standard model, and the Ehrenfest classi- 
fication is used to determine the order of the Electroweak phase transition, analysing 
the effect of the new fields on this physical process. Finally, we show that the new 
fields can greatly increase the stability of Electroweak Z strings and discuss an internal 
space version of torsion gravity sourced by the strings and its effect on the fermion 
masses in addition to that of the Higgs force. Kibble and Lepora |EVGj have studied 
two metrical structures on the Electroweak vacuum manifold. Applying their idea, we 
find a set of scalar metrics, one for each fermion mass. 

The calculation to determine the order of phase transition is carried out for hoth the 
Euclidean and the Lorentzian signatures, and both signatures are studied below in 
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section O Now we explain how the positive features in each lead us to consider them 
each in turn. 

The extension in question was first done in the Euclidean signature |SMVj and Majorana- 
type terms were identified, but the physical meaning of these terms was unclear. In 
|SNMj and [LNGj a modification of the axioms was made, and a way to carry out 
the extension in the Lorentzian signature was realised with the KO-dimension of the 
spectral triple 6 instead of 0. This time, there was a genuine Majorana mass term 
for the uji whose mass was generated via a see-saw-like mechanism, by a new scalar 
field separate from the Higgs. Moreover, the fermion quadrupling problem was also 
solved. Thus a completely vahd physical theory was developed, in particular, as the 
Gelfand-Naimark theorem has no counterpart to non-compact nor noncommutative 
spaces. 

Now we give the reasons why we will nevertheless, also include the Euclidean signature 
case. Firstly, the internal space Einstein's equations of |SMV] led to a fully degenerate 
set of eigenvalues for the fermion mass matrix. However, with the further constraint 
that leptoquarks be defined to be non-dynamical, the correct fermion mass matrix 
turned out to be a solution of the field equations: the solutions had features compatible 
with experiment even in addition to the non-degeneracy of the masses, and for this 
reason it is arguably unlikely that these results came out correctly by coincidence 
|TSMj . Further to this result, we will carry out our calculations in the Euclidean and 
consider the Wick rotation to the Lorentzian signature to take place at the end. The 
reason why these solutions were not included in |SMVj was because the aim there was 
to consider the result of varying all the degrees of freedom in the finite triple Dirac 
operator, in other words no parameters were to be held constant. At first then, it 
appeared unnatural that some fields be non-dynamical. Also, the most general form 
for the Dirac operator was found mathematically instead of the mass matrix being 
added in by hand. Here, we avoid colour symmetry breaking not by adding in a new 
constraint, but by providing a reason why the leptoquarks might be constants. In the 
final section of this paper we explain that the static strings with improved stability 
can source non-propagating internal space torsion, and hence they provide a physical 
reason why the leptoquarks should forfeit their equations of motion. 

Moreover, from a mathematical point of view, a 'noncommutative space' is the gener- 
alisation of a Riemannian manifold to noncommutative geometry rather than that of 
a pseudo-Riemannian manifold, and therefore, in doing calculations in the Euclidean 
signature, nothing of the original spirit of noncommutative geometry is lost; a minimal 
generalistion of the Gelfand-Naimark theorem to an 'almost commutative geometry' 
would maintain local compactness, while the representation of the algebra would still 
act on a Hilbert space, which by definition has a positive definite signature. In addi- 
tion, it has been revealed that there is no Hochschild cycle in the Lorentzian model 
with right-handed neutrinos that can generate the chirality |NOA] . 

The latter is a mathematical argument only, and since the mathematics is not com- 
plete for the i^O-dimension 6 almost-commutative theory (by the arguments above for 
example) anyway, we argue that further study of massive right-handed neutrinos in the 
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Lorentzian signature in the current context should not be ruled out: we might expect 
that a counterpart of the Gelfand-Naimark theorem to Lorentzian almost-commutative 
geometries would lead to equivalent results as those of the ifO-dimension model after 
a Wick rotation. 

In order for this article to make coherent sense, it is necessary to give a minimal 
amount of background information on the articles we cite. Details on the subject of 
the spectral action principle in noncommutative geometry and previous work are kept 
to a minimum. In contrast, we will elaborate a little more on the physics of phase 
transitions in order to motivate the calculations that follow, and also work through 
Dolan and Jackiw's derivation in more detail than is already available in the literature. 

1.1 Background information on noncommutative geometry 

Connes' noncommutative geometry is built by the notion of a spectral triple defined as 
follows, and the spectral action principle is the theory of the standard model coupled 
to gravity whose action is derived by calculating the eigenvalues of the Dirac operator 
that acts on a tensor product space of the spectral triple of space-time and the 
dimensional internal space, finite spectral triple. (Internal space is similar to the extra 
dimensions in Kaluza- Klein theory; it represents the charges on the particles.) As a 
diffeomorphism invariant geometrical theory, the spectral action principle provides an 
important step towards the unification of particle physics with gravity. Moreover, the 
Higgs takes on the geometrical role of connection on internal space, and hence is the 
field strength in the internal space equivalent of gravity. 

A spectral triple {A, TC, V) provides the analogue of a Riemannian spin manifold to 
noncommutative geometry. It consists of 1. a real, involutive, not necessarily commu- 
tative algebra A, 2. A Hilbert space 7i on which the algebra is represented and 3. A 
Dirac operator T) that gives a notion of distance, and from which is built a differential 
algebra. The geometry of any closed Riemannian spin manifold can be fully described 
by a spectral triple and a non-commutative geometry is essentially the same structure 
but with the generalisation that the algebra of coordinates are allowed to be non- 
commuting. In the case of the standard model dimensional internal space spectral 
triple, the Hilbert space represents the fermions, the Dirac operator is the fermion 
mass matrix and the algebra is derived from the symmetries of the action such that 
the latter be the internal automorphisms (lifted to the spinors) of the algebra. 

The Dirac operator encodes the metric, and Schiicker et al (for one generation of 
elementary fermions) calculated Einstein's equations for internal space by minimising 
the Higgs potential - this is the part of the spectral action that depends on the internal 
space Dirac operator only, also called the finite spectral action - with respect to the 
fluctuations of the Dirac operator |OCIj . They found that the standard model fermion 
mass matrix was a solution. In [SMV] the equations were calculated using a different 
method: all the fluctuations of the Dirac operator were varied, including the Yukawa 
couplings and Cabbibo-Kobayashi-Maskawa angles, thus promoting these traditional 
constants to dynamical variables. This was done for two reasons: to apply the idea 
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that in gravity all the degrees of freedom in the metric should be varied, and to attempt 
to remove the need for choice in the pre-fluctuation Dirac operator. To this second 
type of Einstein equation, the solution is that the fermion masses are all the same, but 
it does identify which masses are non-zero and which are zero. 

The standard model spectral triple is a product of two spectral triples, one represent- 
ing (Euclidean) space-time and the other representing the internal space of particle 
charges. The latter has a finite dimensional Hilbert space and a noncommutative 
algebra. The tensor product algebra is 'almost commutative': 

>1 = C~(M)0(H©C©M3(C)) (1) 
The internal Hilbert space is: H = Hl® 'Hr © TY^ © 'H^ where 

Hr = ((C © C) ® (g) C=^) © (C (g) C^) 

and whose basis is labelled by the elementary fcrmions and their antiparticles. The 
symbol c is used to indicate the section represented by the antiparticles. In cither case 
of Hl and TCr, the first direct summand is the quarks and the second, the leptons. 
N stands for the number of generations. For example, the left-handed up and down 
quarks form an isospin doublet and their right-handed counterparts are singlets and 
there arc three colours for quarks and none for leptons. The charges on the particles 
are identified by the faithful representation of the algebra on the Hilbert space. 

Since the fermions can be left or right-handed, H is Z/2 graded and the grading 
operator is given by the chirality operator x = diag(l, —1,1,-1). If % anticommutes 
with V then the triple is called even, otherwise it is called odd. Ti may be written as 
a direct sum of two orthogonal parts represented by particles and antiparticles. If it 
possesses the additional grading operator 7, which commutes with each of T>, x and 
p{a) and whose eigenvalues are +1 for particles and —1 for antiparticles and is called 
(S^-real. An 5'°-real triple excludes Majorana spinors. The standard model excludes 
Majorana particles and is formulated as a real, S'^-real, spectral triple. 

For a real spectral triple, there is an antihnear isometry J. The action of J on as 
given by the composition of charge conjugation and complex conjugation: 




where the bar indicates complex conjugation. 
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J obeys the following relations for KO-dimension 6 and 0: = 1, JD = DJ. For 
i^'O-dimension 0, J commutes with chirality x ^^id in the case of i^O-dimension 6, J 
anticommutes with x- The last two data correspond to the fact that in the Euclidean 
signature the chirality is unaffected by charge conjugation, while in the Lorentzian 
signature it is. 



2 Field content 
2.1 Introduction 

Before we can start the calculations indicated in the introduction, we first have to 
establish exactly what is the field content of the theory both on and off-shell in each 
signature, and underline that the new fields do not conflict with experiment: we check 
that colour symmetry breaking does not occur, and that the additional fields do not 
add any new symmetries to the action. In this part we also show that new fields can 
provide Majorana terms that are non-zero in the physical vacuum. We will treat the 
Lorentzian and the Euclidean signatures in turn. For the latter case, we discuss the 
Wick rotation. The calculations in the two articles referred to in this introductory 
subsection were carried out in the Euclidean signature. 

Paschke, Scheck and Sitarz showed that the new fields they called "leptoquarks" arise 
as those elements of the Dirac operator allowed by the first order condition when the 
5'°-reality condition is omitted. The effect of S'°-reality is to force the blocks labelled 
Z below to in order to disallow Majorana terms as in the standard model. The 
leptoquarks are elements of Z. 



These new fields comprise 6 complex scalars having both the colour and weak isospin 
charges. They are due to the following connections between the quarks and leptons 
in the calculation of the first order condition as a constraint on Dp: This calculation 
has already been demonstrated fully in the papers we refer to, and to avoid needless 
repetition we will not give it here. The generalised Dp without the S'°-reality condition 
is presented in the case of each signature below. 

The generalisation given in |SMVj comprised of omitting S'°-reality, and also including 
the uji fermions. □ Two more new fields appeared for each fermionic generation, 

^As the matrix given here is square, the basis really includes vr, although a row and column of 
zeroes can always be written in for standard model considerations. 

•^To avoid violating Poincare duality the chirality right r neutrino was left out. 




(2) 



where the bar indicates complex conjugation, and the basis is the particles and an- 
tiparticle sections of the finite dimensional Hilbert space. 
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specifically a weak isospin doublet. As clarified below, f l2.3p this extension resulted in 
Majorana-type terms for the neutrino. 



2.2 Lorentzian 



"A Lorentzian version of the non-commutative geometry of the standard model of 
particle physics" was worked out by Barrett [LNG] and similarly by Connes [SNMj in 
2006. In this model the generalisation done in [SMVj by adding the uji and omitting 
S^-reality was carried through to the Lorentzian signature. Since leptoquarks appear 
also in this model, they had to be ruled out to avoid colour symmetry breaking: in 
|GSNj a constraint was applied to the Dirac operator to delete their presence from the 
action. 

From the following axioms [SAP] in dimensions: DpJ = JDp, DpX = ~xDf, 
Dp = Dp, [[DF,a],b"]], [[DF,b"],a]], one arrives at the following complex matrix 
presentation for Dp with Hilbert space basis (\E' = ipL,ipii,ip'i,ip'^), where each ip 
represents the set of standard model fermions with the addition of course of z/^j: 



Df 



f 





M* 








\ 




M 








H 
























H* 


M 





/ 



(3) 



where M = Mq ® I3 © Ml, that is, the Dirac mass matrix for quarks Q and leptons 
L and where H is symmetric. In symbols, H = and H = H* where t and * 
denote transposition and hermitian conjugation respectively. Considering the basis, 
we note that the blocks labelled M and M* are maps between left and right and they 
give Dirac masses by the Higgs mechanism, whereas if is a map from right to right, 
which is indicative of Majorana mass terms. (Note that in the Lorentzian signature, 
J anticommutes with x, and so in this case, x = diag(l, —1, —1, 1)). 

After applying all the commutation rules above involving Dp (plus the new constraint), 
it is seen that for each family of fermions, there is only one non-zero entry which 
appears in the diagonal place in the final row and column (in the convention of [SMVj ) . 
which we label 'p' and it arises due to the connection between chirality right neutrinos. 
When Dp is substituted into the Yukwawa action terms: (\1/, the field p gives 

the Majorana mass term: vr p vr and not surprisingly, from p arises the hermitian 
conjugate term. \i i/r = Q as in the standard model, p is 0. Since p only gives a mass 
to the vr, the vl remains massless. 

Next we check that the field p is non-zero in the physical vacuum. The reason for 
asking this question is that the solutions to the internal space Einstein's equations 
in [SMVj identifies zero and non-zero masses, and if we find that p is zero by this 
method, it would not provide a non-zero mass in the Einstein equations calculated by 
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the method in [OCIj either. So to check that p uuhe a vahd Majorana mass term, 
we make this verification. To do this it is enough to consider the first generation alone. 

As before, to calculate internal space Einstein's equations we vary the finite spectral 
action: @ 



S = Tr(Z}^ - 2Dl) 



(4) 



with respect to the degrees of freedom in Dp. 

Provided there are only quadratic terms or higher order, it doesn't matter whether 
we identify the zero elements of H before or after performing the differentiation. The 
action S is contains no lower order terms than quadratic terms in both H and M. 



S = TT[2{MM*y + {HHy + AMM*HH - 2{HH + 2MM*)] 



(5) 



This was obtained by inputting [3] into [Hand simplifying by using cyclicity of trace and 
the trace of a matrix is the same as the trace of its transpose and finally by writing 
the action as one part plus its hermitian conjugate and taking only the first part. 

Moreover the equations of motion can be calculated just by differentiating H] with 
respect to Dp to get: 



Dr 



(6) 



which tells us that: 



DfD*pDf = Df 
D*pDfD*p = D*p 



(7) 
(8) 



in other words the vacuum is a partial isometry. It has eigenvalues in the set -1, 
and +1 and so the diagonal form of D\ is a diagonal matrix each of whose elements 
are 1 and 0. Therefore we write out the matrix equation D\: 



J^F 



( M*M 



\ HM 





MM* + HH 
M*H 






HM 
M^M 




M*H 



HH + MM* 



\ 



(9) 



^This is obtained by taking one point only in space-time, and setting the gauge fields to zero. 
With Dp the standard model fermion mass matrix, it is equivalent to the usual Higgs potential. 
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It is clear that the field p does not enter the equations involving the quark Dirac mass 
matrix and from the corresponding calculation done in |TSMj it is immediate that 
the vacuum solution of the latter is a partial isometry. We can now work with only 
the first generation lepton basis and so M and H can be taken to be 2 by 2 complex 
matrices. 

--(:^)- "^ol) 

Performing a gauge transformation so that each mass matrix M and H is diagonal, we 
find the vacuum solutions: H and M are orthogonal complement diagonal projections. 
Therefore the following two cases are solutions: 

and 

The first solution given above tells us that the field p (bottom right hand entry of H) 
has a non-zero vacuum solution and z/^j has a Majorana mass and no Dirac mass. In 
the other solution given above, has a Dirac mass and no Majorana mass. Seeing as 
H is a field independent altogether to M, the neutrino mass is unrelated to the Dirac 
masses of the other fermions. 

The pertinent action terms (fTUl below), followed by the resulting equation of motion 
for the field t ( ITTI) are: 

{tif + 2tt (-1 + rf + ss+ pp) + 2fsqt (10) 

i{\t\^ + |r|^ + |sp + \p\^ - 1) + fsp = (11) 

Conclusion: the field content of the theory is that of the original spectral action prin- 
ciple plus the singlet complex scalar p. There are no contradictions to experiment as 
there are no leptoquarks to cause colour symmetry breaking and p generates a Majo- 
rana mass in a similar way to the see-saw mechanism. Also, the introduction of these 
new scalar fields does not change the symmetry group of the action as the terms in 
the singlet complex scalar field p are U{1) invariant. 
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2.3 Euclidean 



All of the calculations above have already been applied in the case of the Euclidean 
signature, so there is little original work to be done for this subsection. Nevertheless 
there are some important clarifications to be made before moving onto the following 
section. 

In |SMVj the was added to the fermions and the S'''-reality condition was omitted. 
As a result, the axioms constraining the finite triple Dirac operator Dp led to the 
appearance of the leptoquarks that had already made their debut in [NGLj plus the two 
new fields j and / that arose from the first order calculation oweing to the connection 
between and uj^, and between z/£ and upt respectively. They form a weak isospin 
doublet. Again we will not demonstrate here the first order condition derivation of 
the new fields because it already well documented. We only state the result of the 
calculation: 



/ M* G \ 

M 

G 

\G* M / 



(12) 



with M = Mq (S> I3 © Ml as before. Notice that G and G* are maps from ip^ to ip'^ 
and i/jr to ip1 respectively indicating Majorana type terms. Unlike H in the previous 
discussion, G is not a symmetric matrix: 



G 



( ° 




















o\ 
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\ y 


V 


h 














I J 



(13) 



where the basis is the set of chirality right antifermions in one family: six quarks and 
two leptons. Considering the basis and the first order condition result, we see that 
the leptoquarks {g, u, x, h, v, y) should appear in the bottom left-hand corner of the 
matrix and that the doublet (j, /) appears in the bottom right. 

By cyclicity of trace and freedom to take the transpose, the finite spectral action H] is 
equal to: 



S = TT{-2{G*G+MM*) + {MM*f + {G*Gf+2{MGG*M*+MM*G'^G+MGMG)) 

(14) 
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We see then that the leptoquarks will appear in the action unlike in the case of the 
Lorentzian signature, but they turned out to be zero in the physical vacuum |SMV] . 
In the symmetric phase vacuum where M = 0, the leptoquarks are non-zero. 

In |TSM] it was demonstrated that if the leptoquarks are non-dynamical, then M in 
the vacuum has many features compatible with the standard model. In particular, 
it has a set of non-degenerate eigenvalues and MM* is diagonal while M*M is not 
diagonal. 

(Since they are constant, the leptoquarks do not bring into effect a spontaneous break- 
ing of colour symmetry). At first it seems unnatural that the leptoquarks be non- 
propagating, but these results seem unlikely to be just coincidental; and in the final 
section of this article, we suggest that the embedded strings produce a non-propagating 
torsion on the internal space and that this torsion can be associated with their stabil- 
ising fields, which are the leptoquarks. 

There are solutions with M ^ and either j or / non-zero so they have the status 
of dynamical fields and they may generate the Majorana type terms eljen, e'j^j*eL 
and i^lli^R, i^jJi^L (plus hermitian conjugates). This shows us that the value of the 
Majorana mass of the must be equal to that of the uji. We also find that this 
mass term for the ul would not occur were it not for the since / is only non-zero 
because is there. These two results seem to be unphysical, but we are working 
in the Euclidean signature, so we have to remember to Wick rotate to the physical 
spinors after the calculations. Since one of the aims of this section is to check for 
physical relevance, we explain that here. 

Unlike the Dirac mass which is a map from left to right, the Majorana mass is a map 
from left to left or right to right: 

mM(*2^L) = mMii^l^'k + 'iPr'^l) (15) 

where the second term is clearly the hermitian conjugate of the first term. It has a 
similar right-handed counterpart. It is clear that the right-handed neutrino ^ r does 
not need to exist for the left-handed neutrino to get a Majorana mass. 

In contrast, a Euclidean Weyl Majorana spinor is given by: l = ipL + V'i- And its 
right-handed counterpart is: r = ipR + ipR- lu this case of the Euclidean signature, 
J commutes with x = diag(l, —1, 1, —1) which means that changing a particle into an 
antiparticle and vice versa does not affect the chirality, and so if we apply the charge 
conjugation operator to "^r or to \I'l we find that the Majorana condition is satisfied. 
The left-handed spinor of the dual basis is: = ^2 + V"! The Euclidean Majo- 
rana neutrino terms given above are equal to the product: ttim'^l'^r plus hermitian 
conjugate: mjv/^/j^L if and "^r are left and right Euclidean neutrino fields. This 
shows that both chiralities must be present for a Euclidean Majorana term. Now let 
us compare the Minkowski and Euclidean terms. 

The Euclidean Majorana terms (which are generated by the map /: 
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rriM'^l'^R = mMiipl^'k + i^l *iPr) 



(16) 
(17) 



if we take the first term from the first row and the second term from the second row and 
perform a Wick rotation then we get [T51 The remaining terms give the right-handed 
counterpart. 

The results from the first order condition calculation show that vl may not acquire 
a Majorana mass if the right-handed neutrino is not present because the existence 
of / 7^ requires that the representation for vr connects with the representation of 
the algebra on the antiparticle leptons. A physical reason for this could be that the 
Minkowskian Majorana (Weyl) spinor as stated above comprises a projection onto 
both spin states and if vr is missing in this Euclidean framework, all the neutrinos 
would be left-handed: vl and its antiparticle z/£ and so the Euclidean neutrino would 
not Wick rotate into a physical Majorana spinor. 

To make a physical interpretation of the internal space Euclidean Majorana terms for 
vl and vr found above, we would make the identification back up to the Euclidean 
space-time spinors (in the same vein as the identification due to fermion doubling) 
and then perform a Wick rotation of the Euclidean spinors to the physical Minkowski 
spinors. (In |LNG] the operator J was redefined for the Lorentzian signature and the 
entire calculation of general D p was done from the beginning in the Lorentzian instead 
of doing the Wick rotation at the end of all the calculations.) 

Conclusion: the field content is the same as that of the original spectral action principle 
plus the leptoquarks and plus the weak isospin doublet (j, /). Hence, the action terms 
in this doublet are SU{2) invariant. The leptoquarks do not cause colour symmetry 
breaking because they are not dynamical. Since j and / have no colour-charge, they 
cannot cause colour symmetry breaking. It is assumed that the leptoquark Majorana- 
type terms have very small couplings. 

3 Quantum and finite temperature corrections in 
(f)^ theory 

In this section we follow the Coleman- Weinberg mechanism [CWj for the calculation 
of quantum corrections in 0^ theory and give details of the Dolan and Jackiw finite 
temperature corrections |DJj in the same theory. The aim is to present a pedagogical 
explanation of these useful procedures for their own sake, while preparing for the 
subsequent section in which we study the Electroweak phase transition in the theory 
that has become referred to as noncommutative geometry with neutrino mixing. 
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3.1 Phase transitions 



Here we provide a general introduction to phase transitions and leave a more detailed 
decription of phase transitions in physics until the opening of the next section. 

A continuous cosmological phase transition occurs when the universe undergoes an 
adiabatic expansion so that the temperature drops below a certain critical temperature 
Tc while a scalar field acts as the order parameter, taking on a non-zero vacuum 
expectation value randomly chosen from the range of possibles that is the vacuum 
manifold. During a phase transition, topological and embedded defects can develop. 
As solitonic solutions to equations of motion in the cosmological context, these defects 
analogous to those found in laboratory experiments also occur in cosmological theories. 
These are discussed in the final section. Phase transitions can occur continuously or 
discontinuously as explained below. 

A general form for the potential F of a real scalar field (f) is: 

where the constants o; to 7 depend on the temperature, pressure and affect the vacuum 
expectation value of the field. The shape of a F — plot is a parabola that narrows 
as (f) increases. In the symmetric phase the expectation value of is zero whereas if a 
becomes negative, (this is the mass term) there is a second solution to the equation of 
motion and gets a non-zero expectation value corresponding to its value in the new 
vacuum, which is the new phase with reduced symmetry. If j3 is non-zero then the 
phase transition occurs suddenly and the transition is called first order, for reasons 
clarified in the next section. If /3 = 0, then the change occurs gradually and this 
is called a second order phase transition. However, for a certain range of physical 
parameters, a theory with no cubic or linear terms can still encounter a first order 
phase transition. 

At small (/), the quadratic term is more important than the quartic term, whereas for 
large 0, the significance is reversed. Due to a symmetry breaking the sign of a changes 
and in the case of a gradual phase transition, the shape of the plot at small 4> changes 
into an upside-down parabola while a new minimum develops at a non-zero value of 0. 
In the case of a sudden transition, the old minimum remains and is separated from the 
new minimum by a potential barrier. The old minimum is called the 'false vacuum' 
and is classically stable due to the barrier, but if the new minimum is at a lower 
potential than the old, then it is energetically favourable for the field to move into the 
new phase and it can quantum tunnel through the barrier to the 'true vacuum'. The 
greater this potential difference, the greater the probability that the phase transition 
happens. 

A first order transition proceeds by nucleation of bubbles of new phase matter inside 
that of the old phase. The bubble walls move outwards and increasing volumes of old 
phase are swallowed up and converted into true vacuum. Colliding bubbles coalesce 
until all of space is filled with true vacuum, and the process is done. 



15 



So, in zero temperature classical field theory, phase transitions do not take place 
because there is no spontaneous symmetry breaking mechanism to drive them. The 
mass term does not change sign so the vacuum state of the field remains zero-valued. 
In practice however, phase transitions as outlined above, do occur due to spontaneous 
symmetry breaking. This means that the classical picture is not complete and it is 
necessary to perform corrections to the potential by considering the zero temperature 
quantum effects. Physically, quantum fluctuations of the vacuum may alter the form of 
the potential so that the field can develop a non-zero vacuum expectation value. As a 
result the vacuum loses a symmetry and the matter obeys a smaller set of conservation 
laws. At finite temperature, (T ^ 0) phase transitions can be brought about by 
temperature fluctuations. The two cases will be considered here in turn. 



3.2 The Coleman- Weinberg Mechanism 



In this subsection we will outline the calculation of the one-loop correction at zero 
temperature to the self-interaction potential of a single, real scalar field (p and show 
that it leads to symmetry breaking. The need for the modification is due to quantum 
fluctuations: (f) ^ (pd + rj. 

Spontaneous symmetry breaking resulting from quantum fluctuations can induce a flrst 
order phase transition if it is implemented by quantum tunnelling - not a continuous 
process. The kinetic (derivative) terms of the Lagrangian are unaffected by these 
corrections. The alteration is to the interaction terms and the result is called the 
effective potential, Ve// = Vciassicai{4>) +0{h) Here we will follow the work of Coleman 
and Weinberg who computed these corrections for a real scalar held (p^ theory. The 
process of symmetry breaking studied here is called the Coleman- Weinberg mechanism. 

The theory we will consider throughout this part is that of a real scalar field. The 
potential is: 



V = im0^ + ^<P' 

Whose value is the free energy density with the expectation value < > of the field. 
< > is the 'order parameter' which will be defined in the next chapter. To calculate 
the quantum fluctuations we perform the substitution — > 4>ci+fl into the action. The 
generator of Green functions is Z = e^'^'-^' Where W[j] is the generator of connected 
Green functions. Completing the square simplifles the subsequent calculation. The 
result is: 



lw\j] = ^ (^SM + J j<Pd^ - ^ IV ln(n + m' + ^<f>l) + 0{h) 
=> W\j] = S[U + j J(i>o + Tr ln(n + + + 0{h') 
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Now we perform the Legendre transformation: (which is famihar from switching from 
Helmholtz to Gibbs free energy.) Seff = W[j] — J j(j)ci and we obtain: 

Seff = Sci + '-^ Tr ln(n + + ^) + 0{h^) (18) 

Since (p^i is a constant, the derivative terms vanish and as the action is given by 
S = J Ld'^x, the effective potential K// = ~yT^^f} where VT is the volume of space- 
time. The final term is of course negligibly small and we will ignore it from now on. 
The infinite series of Feynman diagrams are obtained term by term from the expansion 
of the logarithm. 

It is impossible to compute the full effective potential because that would require a 
summation of an infinite number of Feynman diagrams. Instead we can either take a 
low order of perturbation theory and compute as many of those diagrams as possible 
or we can select one or two typical diagrams from each of many orders. To lowest 
order, (called the tree approximation) the Feynman diagram is a cross: it has four 
legs and a vertex. It indicates the original term ^4>^- The four legs represent the four 
powers of and the vertex represents the coupling constant. To first order, (one loop) 
there is an infinite sequence of polygon graphs begun below. Their sum is called the 
one-loop radiation correction. 

The number of external legs represents the number of powers which shows the 
number of particles taking place in the interaction. For each vertex we write A and 
integrate over the four-momentum. We include the | to account for bose statistics; 
this is to prevent us from overcounting equivalent diagrams. Each internal line carries 
a massless propagator p^^. A factor of ^ is also to prevent overcounting: reflection 
or rotation of the polygons leads to the generation of equivalent diagrams, not new 
contractions in the Wick expansion. 

The resulting interaction terms of the effective potential are as follows. The second 
term is the coupling constant counter term with coefficient C . The ubiquitous h of 
quantum theory also appears with the i as a prefactor but is not shown here because 
we are working in natural units where h = c=l. 




The final term in [19] (we will label as 'Vi') looks like the formula for the expansion of 
a natural logarithm ln(l+x) for small x. Transforming to Euclidean space, ik^ — ^ k^, 
we rewrite it as: 

T. I f d^k ^ . \(P\ 
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By dimensional analysis (or power counting) we can see that it is divergent. Coleman 
and Weinberg performed a cut-off regularisation at k"^ = which gave: 



647r2^ 2567r2 ^ V2AV 2 ' 

It then remained to renormalise the mass. (The resulting expression depends on 
the choice of renormalisation conditions. The choice does not alter the physical the- 
ory.) When this was performed, it became necessary to shift the fields which in turn 
prompted a redefinition of the coupling constants. This led to the determination of A 
in terms of the coupling strength e. The final result was: 

where a is the renormalisation scale. 

With the inclusion of the radiative correction order h term Vi into the potential it 
becomes V^jf = Kiassicai(0) + 0{h) This quantum correction drives a breaking of the 
symmetry of the Lagrangian. Without this term the the solution to the equation 
^ = is 101 = but the solution to ^^^^ = gives a non-zero expectation value of 
the field. 

Although the steps outlined above were carried out in the case of a single, real scalar 
field, the methodology can be applied to other theories. For a non-Abelian gauge the- 
ory with scalar fields, gauge fields and fermion fields, the Coleman- Weinberg correction 
with the 'zero-momentum' renormalisation conditions is as follows: 



V^i(0) = -— [ Tr ( fiHn^ ) +3Tr (MMn — ) - 4Tr ( m^ln— ) ] 

where/i, M and m are the scalar, vector (gauge fields) and spinor (fermion) masses 
respectively. 



3.3 Finite temperature corrections 

At finite temperature, temperature fluctuations are clearly far more significant than 
quantum fluctuations. Corrections to the potential must be made to account for the 
interactions between particles in an ideal gas at finite temperature. Dolan and Jackiw 
calculated the finite temperature corrections to the potential. 

With the aid of [R] we will reproduce their derivation here in detail. 

Firstly we recall the result of the first order summation of the Feynman diagrams from 

m 
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Vi = lJ ^ln(fc^ + f/"(0)) (21) 

which is divergent. Note that U" = Ai^, that is, the coefficient in the quadratic term. 
{h is set to 1). 

By switching to Euchdean space and integrating out the zeroth component of the 
momentum (this step is shown below), we obtain: 

^[-a; + Tln(l-exp(-/?^))] (22) 
where (3 is the inverse temperature. 

In the zero temperature hmit studied above, the second term clearly vanished (T = 
0) and the divergent ffist term was analysed. The second term which is due to a 
mixture of quantum and temperature fluctuations is not divergent. In this case of 
high temperatures, we omit the first term which represents the zero point energy. At 
zero temperature, quantum fluctuations are insignificant compared to variations in 
temperature so that changes in the ffist term are much smaller than changes in the 
second, although some authors have considered both terms simultaneously. 

The steps to be taken to derive equation [22] above from [21] are as follows: 

We impose the following boundary condition upon the Euclidean time coordinate: 
'Time period'=2^ where T =temperature and h is set to 1. 

In this step, we have not replaced physical time with inverse temperature, and there 
is no physical interpretation of the relation. After all, the Euclidean time coordinate 
is imaginary, having no physical significance. The quantities to calculate in thermal 
quantum field theory are thermal expectation values of observables, < A >. These 
are normally independent of time and so the assignment of time dependence does not 
affect the physics. This provides a freedom of choice that can be exploited by treating 
the density matrix as a time evolution operator. This is explained below. 
The time-ordered thermal Green functions are defined as: 

iGi''\h,...,t^;T) = ^Tr{e-^"T[Mri)-MrN)]} 

The thermal expectation values differ from those calculated in zero-temperature quan- 
tum field theory as the time ordered product involves the density matrix, e~^^ and the 
trace (Tr) over all the possible states in the system. This brings the thermal average 
of statistical physics into the quantum theory. The density matrix looks like the time 
evolution operator in Quantum Mechanics. If we make the replacement: 
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"Ti—if} 

e'^" ^ exp f - / drH 



(3 defines the periodic boundary conditions on the imaginary time coordinate. Let 
Tj = then < r < /5. For fermions, the corresponding boundary conditions are 
anti-periodic. 

ko = Nfku = /i V — = 27r V27rT 

' time period ' 

n n 

One of the factors of 27r cancels with that in the invariant measure of Vi when we 
integrate out the kc,: 

n=+oo 



dko — > 2nT ^ 



With this transformation the finite temperature correction becomes: 

d^k 



where a;|„ = {27mTf and kf + ^ E"^ 
Let 



n 

dE ^W + ^L 



+ (27rn/3-i) V E E^ + {27mp-^y 



2 4/3^ §^ 2 2/3 /-TT 71 ^ 

^ + E^B^ ' , ^2 = ^ + + ^COth( — 



+ e"'^^ _ /3e-^(e2^ + 1) _ /^(e^"' + 2-1) 
- e-^ e-^(e2^ - 1) ~ e^^ - 1 



2 
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since x = ^ 



Let us restate this result more clearly: 



dE V2 - 1 

/;jl|dE=„(E)=/2/j(i+^)dE 

= 2/3|i£; + /3-'ln(l-e-'«')) 

which is the result we were looking for. The last step can easily be checked by carrying 
it out backwards. That is to say differentiating the last expression gives the integrand 
in the penultimate. 

The final procedure is to write [22] in a form in which we can easily read how these finite 
temperature corrections can effect a symmetry breaking. If we assume that A4. 
we can perform an expansion in powers of These steps are set out below. 

First we introduce the function I{y) = dx x"^ ln[l — e^^'^^'^^^^] where y = (3Ai and 



so Vi 



d^k — > dx, so to ensure we do not make a dimensional error we include the x"^. 
Taylor expanding: 

I{y)=I{0)+yI'{0) + ^L_^ + O{y') 

Evaluating the terms: 

m 



Iq = dx x^ ln(l - e"^) = - rT^ / dx x^ e "~ 

n=l J 



45 



dl{y) 



dy 



dx x\\ - e-^^ +y ^~'Y^-2y{x^ + y^-^e-^^^^ ^ 
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evaluating at = . . . 

oo -| poo 

dxx\l- e-'^y^-e-'' = / rfxxe-^(l-e-^)-^ 
X Jo 

2 ^ 2 Jq ^ ' VI 

we substitute these back into V\ = j:^^ to obtain the temperature dependent contri- 
bution to the free energy: 

Vi = -—T^ + -—T^ + 0(T) (23) 
90 24 ^ ^ ^ ^ 

with the replacement y — > 

The 0{T) term is neglected in the high temperature limit. 
The final effective potential becomes: 

VefM T) = -^0^ + - ^T^ + (24) 

The physical interpretation of the last result now follows. 
With _ = - m2 

the last equation YIM becomes: 

^1 = + — - = + 02 

^ 90 24 ^ ^ ^ 90 8 ^ 24 

m202 AT202 



(25) 



-|-terms not dependent upon 0. 

The field varies smoothly as the temperature changes. A temperature change from 
one side of the critical temperature to the other would induce a second order phase 
transition. For example, the equation holds that if the temperature is increased up to 
the critical temperature, the field decreases continuously to zero. 

In general the 02 term yields the mass squared. Inspection of this coefficient yields 
crucial information about how a phase transition is implemented by a temperature 
change. 

- = - ^'W = (26) 

Four points to consider about the above expression are: 
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• When JH^ < the field obtains a vacuum expectation value and the symmetry 
is broken; T < Tc 

• When j\4^ > the symmetry is maintained and = 0; T > Tc 

• When the temperature T increases from below T^, M.^ becomes =0 and then it 
becomes positive and symmetry is restored. 

• The critical temperature is defined by = 0. 

It is evident from the above equation that the approximation T ^ we used in 
the expansion carried out in the previous subsection is valid provided the coupling 
constant A is small. 

4 Analysis of the effective potential: determining 
the order of the phase transition in two theories 

In this section we will demonstrate how to determine the order of a phase transition 
from the effective potential in finite temperature field theory. In the previous section, 
the finite temperature corrections were demonstrated only in 0^-theory, but equation 
121] can be applied to any theory just by taking the second differential of the potential 
in that theory: U" . It follows that we may write down the effective potential for 
noncommutative geometry with neutrino mixing, and subsequently analyse it in the 
Ehrenfest classification (see below) in order to predict whether the Electroweak phase 
transition was to be first or second order. We also assume that after a Wick rotation 
back to the Lorentzian, we can still add on the correction terms in the same way 
where the meaning of T becomes physical temperature, relinquishing its connection 
with imaginary time. 

The 0^ theory is used like a control experiment as we already know that the answer is 
that it be a second order phase transition because does not have a linear or cubic term. 
And these derivations are carried out side by side for a clear and direct comparison. 
The purpose of this section is to carry out a corresponding calculation in the case of 
noncommutative geometry with neutrino mixing to ascertain whether the new fields 
in the theory change this. The Electroweak phase transition is first order for a certain 
range of physical parameters, but we ask whether the new fields ensure that it be first 
or second order for a greater range of parameters. Of course, the new interactions are 
quadratic in the Higgs, so if the leptoquarks are non-dynamical, we expect the phase 
transition to be second order. However, since as dynamical fields, they interact with 
the Higgs in the vacuum, we should work through the calculation to find out whether 
they have any effect on the order of the phase transition. 

First we explain a little about the relevant physics. 
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The Ehrenfest classification (1933) of phase transitions states that if the nth deriva- 
tive of the chemical potential /i is the lowest order discontinuous derivative then the 
transition is nth order. We will discuss ^ rather than ^ since the Gibbs free energy 
G and chemical potential are related in the following simple way: 



For a first order phase transition is discontinuous. The heat capacity at constant 
pressure Cp is given by the derivative of enthalpy H with respect to temperature. This 
is also infinite since the enthalpy changes a finite amount for an infinitesimal change in 
the temperature. For example, water heated to its boiling point remains at the same 
temperature as heat continues to be applied to it since the energy is being taken up 
by the bond breaking rather than to increase the kinetic energy of the particles. 

In general, if there is a non-zero latent heat released during the phase change then the 
transition is first order. When heating a liquid until it evaporates its pressure maintains 
a constant value p < Pc where pc is the pressure at the critical point. When bubbles of 
the new phase form, the pressure jumps to its new value. Ap has a simple and direct 
relationship with the change in free energy. If the evaporation experiment is repeated 
but this time the pressure is held at constant pressure p = Pc the latent heat becomes 
smaller and in the limit, no heating is required to convert the liquid into vapour. At 
the critical point, the latent heat of vaprization is zero and as a consequence, the first 
derivative of the Helmholtz free energy with respect to temperature is continuous at 
that point and so we have a second order phase transition. 

In the case of a second order phase transition the first derivative of the free energy 
with respect to temperature is continuous but the second order derivative is not. This 
time the heat capacity at constant pressure changes discontinuously but is not infinite. 
An example of a second order phase transition can be given as a semiconductor whose 
temperature increases beyond the critical temperature and it becomes an insulator, or 
the recrystallisation of a metal by annealing. 

Below we will use the Ehrenfest classification to analyse the effective potential in 
finite temperature field theory and compare the forms of the potentials for first and 
second order phase transitions. 

4.1 Analysis of the effective potentials 

The following cases are considered in parallel: 

(A) : ordinary 0^ theory 

(B) : noncommutative geometry with neutrino mixing in the Lorentzian signature. 
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(C) : noncommutative geometry with neutrino mixing derived from an irreducible (in 
the sense of |SMV] ) spectral triple and with colour charge neglected, in the Euclidean 
signature. 

(C2): the same as C but with constant leptoquarks. 

(D) : the same as C but in the theory with the full fermion basis to include their colour 
charge. 

{D2): As [D) but with constant leptoquarks. 

The reason that the cases {B), {D) and (-D2) are each studied is explained in the 
introduction and the notation is defined in [21 Cases (C) and (C2) are also worked 
through in order to make the comparison to analyse the effect if any, that the existence 
of colour charge may have on the derivation. Case [A) is included to act as a control 
experiment with which to compare the stated theories. For some formulae and some 
steps, (C2) and (-D2) are the same as (C) and (D) respectively, so we do not always 
write them separately. 

In writing down the potentials V(^) to V(£)2), [fl we will study only one row from each 
set of matrix equations that are the equations of motion (integrated) derived from the 
finite spectral actions presented in [2] (for example [TT]) because for the purposes of the 
calculations to follow, it is enough to consider one row at a time. Not all the details as 
to the derivation of the potentials are given here because this information has already 
been presented and analysed in |SMVj . The solutions pertaining to the symmetric 
phase extremum of the potential and that of the new phase, or physical vacuum are 
written down in the tables below. For (B) we will mainly focus on the solution where 
p is zero in the physical vacuum. The field labels, t, a and so on have all been defined 
in [21 The equations of motion in the Lorentzian case were given earlier in [21 whereas 
those for the Euclidean case were calculated in articles cited in that section. Let / 
denote the symmetric old phase and // the new phase with less symmetry. 

V^B) = \t\^ + 2|t|'(|r|2 + + - 1) + 2fsqt 
V^c) = + 2|ap(-l + \b\^ + |cp + 1^1^ + \h\'^) + 2cbda 



V^D) = 3|a|^ + 6|a|2(|6|2 + \c\'^ + \g\'^ + \u\'^ + + \h\'^ + \v\'^ + \y\'^ - 3) + Gcbda 
^We need only to include the terms in a, t. For example, linear terms in a are excluded. 
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Symmetric phase / | 


(A) 


|0P = O 1 


(B) 






(C) 




= 0, |5p+|/ip=i 1 


{C2) 


|a|2 = |6|2 = o, c/, /i constant, a(|ap + + |c,|2 + |/i|2 _ i) = | 




|ap = |6|2 = 0, \g\' + \ 


up + + Zip + + = 3 1 


m 


|a|2 = |6|2 = 0, 


gi, li, X, h, V, y constant 





New phase // | 


(A) 


1012 = M! 1 


(B) 




|S|2 = 1, |p|2 = 1 


iC) 




|a|2+|6|2 = l, c, = /i = 1 


(C2) 




(D) 


|a|2 + |6|2 = l, 


g=u=x=h=v=y=0 


m 


|ap + |6p + |c,|2 + |-u| 


^ + kl 


^ + |/iP + |T;|2 + |y|2 = 3 a^O, |a|Vl | 
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4.1.1 Finite temperature corrections 



Now we test for what type of phase transition it may have been that happened as the 
Electroweak epoch progressed into the electromagnetic epoch of today. In the way 
already explained, we write down the finite temperature corrections for each case: 



Veff (A) = -r r + — + 

eff (A) ^V^ 2 ^ 90 24 



2 KA2 rT.2 

Veff (C) = + 2|a|^(-l + 16^ + \c\' + \g\' + \h\') + 2cbda - -T' + 



Veff p) = 3|a|^ + 6\a\\\b\' + \c\') + 2\a\\\g\' + \u\' + \x\' + \h\' + \v\' + \y\' - 3) 
where 



2 _ Vl^A) 



and so on is given by: 



MIb) = 4|t|' + 2|r|2 + 2|s|2 + 2\p\^ - 2 



Mj^) = 4|a|2 + 2|6|2 + 2|c|2 + 2\g\^ + 2\h\'' - 2 



Mfn) = 6(|a|2 + + |c|2) + 2(|y|2 + + + \h\'' + + _ 3) 
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4.1.2 Critical temperature 



Substituting for in K//, collecting together the terms quadratic in and omitting 
terms in T but not in (f), t, or a respectively to find the mass terms, while working in 
mass eigenstates, that is, r — s — b — c — 0: 



47^2 

Veff iB) = \t + 2\t\\-l + \p\') + 
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47^2 

Veff (C) ^\a\' + 2\a\\-l + \g\' + \h\') + — 



6T^ 

Veff (D) = 3|a|^ + 2|a|^(-3 + \g\^ + \u\^ + \x\^ + \h\'^ + \v\'^ + \y\^) + — 



with masses: 



l-r2 



-2^(C) = -2(4-4M -4|/ir- y) 



-^M(% = -^(6 - 4|^|2 - A\u\^ - A\x\^ - A\h\^ - A\v\^ - 4|y|2 - ^) 

As the temperature increases, one can see that the mass term goes positive, where the 
only vacuum solution for the field is 0. The critical event, where the mass squared is 0, 
is when the universe is at the critical temperature, Tc- For constant G, the contribution 
to the critical temperature from the new fields is to decrease it, as — jg'p — |itp — — 
— |f p — |yp + 3 takes a values dependent on the square of the mass of the up 
particle. However, it is difficult at this stage to state anything definitive about actual 
numbers. 
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4.1.3 The Ehrenfest classification 



Since we are making a comparison, we may omit any terms that occur identically on 
both sides of the critical event. Here, let us substitute with = to write down the 
potential (using 1^ 1 . 1 1 above) during the symmetric phase: 



Vi, (A) = Vi, (B) = Vi, (C) = ^i, (D) = 



Vi, iC2) = -^\au? = ^{\g? + H-l) 

Vi, = -^3|a,,p = ^ma? + + \x\' + \h\' + \v\' + \y\') - 6) 

At Vii the new phase, we substitute for the non-zero vacuum expectation value of the 
field: 

The formula + \h\^ is interpreted as the Higgs doublet term |0p = + |02p and 
as the true vacuum the Higgs chooses a direction at random we can take 02 = 0, so 
we choose 6 = 0. Similarly, let s = 0. 

V - 
Vn, (A) - 

Vii, (B) = Vn, (C) = Y' = X 

I |2 I I |2 

^11, (C2) = "^1*^1 + (C2) = 

Vn, iD2) = \a\\3\a\^ - 3\a\^) + ^(6|a|') + Vj, (C2) = Y^(3|a|') 
So, Vii^ (C2) = —Vi^ (C2) and Vn^ (d2) = —Vj^ (d2)- 

Next we compare the first differential of the two potentials with respect to the tem- 
perature evaluated at the critical temperature: 

dVi, (A) _ dVi^ (B) _ dVi, (c) _ dVj^ (g) _ 
dT ~ dT ~ dT ~ dT ~ 



At T = T, the modulus of and is equal to that of and 

respectively. 
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At the critical temperature: 



dV, 



dT 



dV, 



II, (c) 
dT 



T 



7^0 



dVi 



II, (A) 



dVn dM^ 



dT 



dM^' dT 



M^T 
4 



0; 



dV, 



II, (D) 

dT 



Tc 



2*° 



The results show that the potentials for (B), (C) and (D) have a discontinuous first 
derivative at T = and we conclude that the phase transition is first order. Whereas 
for (A), (C2) and (-D2), the first derivatives are continuous at the critical event, which 
means that the phase transition cannot be first order. For a second order phase 
transition, the second order derivative of the potential with respect to the temperature 
is discontinuous at the critical temperature. ^^t^^^^ 'dT^^^^ have the same 

modulus but opposite signs to ^'J^^^^^'' and respectively. We expect Vn to 

decrease as the temperature decreases and V/ to increase as the temperature decreases, 
and since T decreases as time progresses, while the expectation value for the field 
increases, for a second order phase transition we expect the second order differential 
with resect to temperature to behave in the same way as the second order differential 
with respect to the field, which is, for the original minimum to evolve into a maximum, 
and the new minimum to become progressively deeper until the phase transition is 
complete. 

In the 0^-theory, the second differential of the potential with respect to temperature 
is discontinuous at T = T^: 



d^V, 



II, (A) 

dT^ 



m 



ATM ^ 



d^Vi, (A) 
dT^ 



4.1.4 Type of extremum in the symmetric phase 

In this part we will establish whether the turning point at = is a maximum, a 
minimum or a saddle point for each of the cases. To do this, we simply differentiate 
V(^) to V(£)) given in part 14. II twice to see how the gradient changes at that point. 



jj,2 U=o = TT 



< 
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dtdt 

d.^Vf 



(^)iM,=o = 2(|pr-l) 



dada 



^^|m.=o, g=i ^ 2{\g\' + \h\' - 1) ^ 



d'ViC2) 
dada 



\Mq: 



=0, constant G = '^{\g\ + 1^1 " 1) < 



The last line is not equal to zero because a = at /, and so the a doesn't cancel in 
the equation of motion: a(|ap + [(^p + — 1) = and since g and h do not have 
their own equations of motion, \g\'^ + \h\'^ — 1 ^ Q . 



dada 
d?Vi 



^|mc,=o, g=i = ^{\g? + \u\^ + + H + l^r + \y? - 1) = 



dada 



Mq=0, constant G = 6(|^| + \u\ + |x| + \h\ + l^l + |y| - 1) < 



In the vacuum where t = and |pp = 1 in the new phase, and p = while |tp = 1 
in the old phase, we get the same potentials as in [B] except that the expressions 
for Vi and Vn are traded. And the two procedures given above can be used to show 
that the phase transition is first order. We can say something new about the critical 
temperature. Recalling that p is the generating field of the Majorana mass, it is 
expected to take a large value, and hence its effect on the critical temperature is to 
reduce it significantly. This is further confirmation that the phase transition is weakly 
first order. 

These results show that the extremum at / is a maximum in cases (A), (C2) and (-D2) 
whereas it is a saddle point in the cases {B), (C) and {D). The former confirms that 
the greater the values of the constant leptoquarks, the greater the range of physical 
parameters for which the phase transition is second order. The latter confirms that 
in the corresponding cases, the phase transition is very weakly (in fact only just) first 
order. Comparing (C) with (D), we find that when colour charge is included, as there 
are even more new fields, the critical temperature is reduced even more, and that the 
maximum at / is even steeper. It is [B] and (C) that are referred to in the abstract 
as these cases are already better understood and recognised than (C2) and (-D2). 



4.1.5 Summary 

The calculations in this section led to the following predictions: 
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(;A^-theory 


Euclidean 


Lorentzian 


Euchdean, 

(t COllStfUlt 


maximum 


saddle point 


saddle point 


maximum 


2nd order 
phase transition 


just 1st order 
phase transition 


just 1st order 
phase transition 


2nd order 
phase transition 



It should be pointed out once again that the scope of these results is relative to the 
range of physical parameters employed. A more extensive study involving the space of 
physical parameters might lead to a more quantative estimate as to what extent the 
new fields can increase or decrease the likelihood that the Electroweak phase transition 
be first order. As mentioned, for a certain range of parameters it is already established 
in the hterature that the Electroweak phase transition was probably first order. The 
aim of our study was to determine whether the effect of the new field interactions 
would be to make it more or less likely that this prediction be correct. 

5 Embedded defects and torsion 

In this section we show that the additional fields in the extension of the standard model 
in noncommutative geometry significantly improve the stability of the Electroweak Z 
string. Using an analogy from condensed matter physics experiments, we suggest that 
the Electroweak Z string solution in this extension is a source of torsion in the internal 
space dimension. Since the string is a static solution, this torsion is not expected to 
be dynamical, and this scenario is given as an explanation as to why the leptoquarks 
should not have equations of motion. The motivation for having constant leptoquarks 
was explained in the introduction. The reasons for introducing torsion in this way is 
that (a) the leptoquarks satisfy the formulation of torsion suggested below, which is 
constructed in analogy to that of Cartan, and (b) because of an argument given below 
in the context of internal space gravity, also called the Higgs force. Finally, we discuss 
the relationship between the energetics of the string and internal space curvature and 
the metric on the vacuum manifold and find that there is a discrete set of these metrics, 
and hence string mode, for each fermion mass. 

5.1 Introduction to topological and embedded defects 
5.1.1 What £ire topological and embedded defects? 

A soliton is a non-dissipating solitary wave. In other words, it "keeps itself together" 
like a pulse sent down a fibre optic cable whose energy remains localised due to the 
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refractive index of the glassj^ Sohtons are particle-hke solutions of the equations of 
motion. Topological defects are solitons with a non-trivial topology to which they owe 
their stability. They consist of localized trappings of old phase matter persisting on 
into the new phase, so they are a relic of symmetric phase matter surrounded by a 
background of true vacuum. 

In a theory that does not allow for such non-trivial topologies (that is, where the 
vacuum manifold is simply connected) it is possible for scalar or fermion fields to 
interact with the soliton, called an embedded defect, and thus increase its stabilityJil 
For example, in the Weinberg-Salam theory of Electroweak unification, embedded 
string and monopole solutions exist but they are unstable to small perturbations. 
This means that they are classically stable objects but they must decay quantum 
mechanically very soon after the phase transition in which they are formed. If this 
were all there were to say, they would be mere mathematical curiosities uninteresting 
to physicists. However, many authors for example Watkins and Vachaspati [BESJ 
have shown that by extending the standard model to include extra scalar fields, the 
stability of the Electroweak Z string may be improved. Now that there is a solid reason 
to make this type of extension, as we have discussed, the extra fields appear naturally 
in noncommutative geometry with neutrino mixing, it is appropriate to consider their 
ideas and to apply them in our context. 



5.1.2 Why do they form? 

Consider the universe in the new phase. The order parameter cannot be isotropic 
across distances greater than that limited by causality: the field takes on different 
directions in different regions at random and it cannot be aligned at points separated 
by distances greater than light could have travelled since the beginning of the phase 
transition. Going back to the analogy with the ferromagnet, domains of a given spin 
orientation form after the symmetry is broken, and defects clearly occur at locations 
where the direction changes (if the spins change from up to down, at some point 
in between, they go to 0) and so topological defects are expected to proliferate near 
domain boundaries. The size of these boundaries is roughly uniform, and is determined 
by the speed of sound instead of the speed of light. The more quickly the phase 
transition takes place, the greater the number of defects form, as there is less time for 
communication between regions. This is why we observe more defects in a quenched 
specimen than an annealed onel^ 

Due to this experience in the laboratory, one realises that if topological defects can 
form, they will form and this fact should be true in cosmology. (Of course, the analogy 
between topological defects formed in condensed matter physics experiments and their 
cosmological counterparts is not exact due to differences related to friction, relativistic 

^There is a sculpture in Lisbon that emits solitonic water waves 

^there is a type of embedded string called a semilocal string that is completely stable, but it does 
not occur in Weinberg-Salam theory nor in the extension we are studying. 

^quenching can improve hardness, which is why the steel is quenched to make a Samurai sword. 
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effects, gravity, causality governed by ligiit speed instead of speed of sound, knowledge 
of initial conditions and so on, but it can still be a useful one to make.) And, many 
grand unified theories predict that the solitonic solutions to the equations of motion 
in those theories must have formed in nature. Topological defects in cosmology can 
trap enormous amounts of energy because their cores contain old phase matter. Recall 
that in the symmetric phase the field is at the top of the potential hill and therefore 
the energy per unit volume will be much greater than that in the new phase. So, if 
they would not dominate the universe, they should at least have been detected by 
telescopes and show their signature in the cosmic microwave background data. As 
mentioned in the introduction, their absence is explained mainly by the fact that as 
the universe expands, distances between points increases and the defects could be 
already so far spread out that not even one of them occurs within our Hubble radius. 
Also, topological defects with opposite topological charge can annihilate cachothcr, 
by unwinding one another, and some cosmic strings can "radiate away" by emitting 
gravitational waves. 

5.1.3 How is torsion sourced by line defects? 

Cosmic strings and indeed embedded strings are line defects. They are analogous to 
the line dislocations and line disinclinations that are associated to broken translational 
and rotational symmetry respectively. The latter is a line source of curvature which 
can be pictured as a removal or addition of a wedge. The former is an edge or screw 
dislocation which are a source of torsion. Consider a particle in the material that 
transcribes a path around the dislocation. The result of this parallel transport is 
characterised by torsion whereas in the case of a disinclination, the result of a parallel 
transport around the defect is characterised by curvature. Therefore, an Electroweak 
string as a static configuration can be a source of non-propagating torsion. 

5.2 Electroweak Z string solution 

Replacing the Electroweak gauge fields, and recalling the spectral Lagrangian (from 
Connes' spectral action formula) with space-time curvature R still set to 0: 



and the equation of motion for the scalar field 0*: 



{d^^ieA^){d^-ieA^)(t>--(t>{\(t>? 




or, using the Lorentz gauge d^A^ — 0, 
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□0 + 0(eM2--(|0p-r^2)) = O 



where □ is the d'Alembertian operator. 

The equation of motion for the degrees of freedom in Dp-, for example t, are found 
similarly, but using the Euclidean Euler-Lagrange equations and hence there is a sign 
difference compared with the above equations of motion. 

- AUt + t{e^A^ + \t\^ + + - 1) = (27) 

where the Higgs potential —1(101^ — rfY has been replaced by + |sp + |pp — 1)^. 
Normally t] is set to 1. When p = then t and s are just the usual Higgs doublet and 
the equation of motion is no different from 15.21 except for the signature change. 

In the Weinberg-Salam theory, the Ansatz for the Electroweak Z string in polar coor- 
dinates is 0(r, 6) = f{r)e^^(t)o where F is a (broken generator of the symmetry group 
of the action in whose Lie algebra the gauge field A takes values. (It has the same 
configuration as the Nielson-Olesen string, but it is not the same string as the latter is 
a topological defect in theories where the vacuum manifold is not simply connected). 
(The gauge field A is part of the Ansatz but since its equation of motion is not given 
above, it is not given here). 0. We can set up boundary conditions so that /(r) — > 1 
as r — > oo and we understand the statement of [EVGj that the boundary conditions of 
the Electroweak strings define geodesic paths on the vacuum manifold because the lat- 
ter are given by 7y = e^* where t is some real parameter. This discussion is developed 
in the final subsection. 

5.3 Stability of Electroweak strings 

Many authors have considered adding scalar fields to extend the Weinberg-Salam the- 
ory and carried out an analysis to determine whether their presence can increase the 
stability of the Electroweak Z string. There are two strings in the Electroweak model, 
the W and the Z. The former is unstable for all parameters while Z is stable for a 
range of values, which in calculations done so far only include unphysical values. We 
will only discuss the Z string. In [EES] . Watkins and Vachaspati add a singlet scalar 
field X with global U{1) charge and describe a "bag" mechanism as follows. The Higgs 
gives X ^ mass by the usual Higgs mechanism, but the back-reaction of x on the Higgs 
is to try to prevent it from acquiring its vev. In other words, the scalar would rather 
live in a region where Higgs vanishes since the mass of x is zero wherever the Higgs 
is zero. The scalar likes to accumulate on the string and tends to maintain the string 
configuration. The string is a bag in which the scalar prefers to sit in and, hence, hold 
together. 

^choice made in |EVG| to minimise the magnetic energy 

^^Kibble and Lepora impose a condition so that there is no flux in the r directions 
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In the current theory, the extra fields are not singlet scalars with a global U{1) charge, 
but as described in [21 p is a complex scalar singlet with local U{1) charge, and (j, /) is a 
complex scalar doublet with weak isospin charge. The action terms in the leptoquarks 
(case (D) in the notation of the last section) have weak isospin and colour charge. 

Recall that H and G are non-zero in the false vacuum where M = 0, and are in the 
true vacuum where M ^ 0. Also recall the equations of motion for H and M calculated 
in [21 Without giving too many details already published, we recall that in [TSMj it was 
shown that G and M are mutually orthogonal complementary diagonal projections. 
(The vacuum solutions in either signature are that Dp is a partial isometry and if we 
assume that M, G are diagonalisable, - otherwise they would not have any physical 
meaning as a mass matrix - and working in mass eigenstates, then we can take the 
solutions to be diagonal projections). This means that the equations of motion for G 
and M or H and M precisely embody the bag phenomenon. 

Watkins and Vachaspati describe both a condensate and a bound state. A bound state 
consisting of the string and the p particles can form as the bag, whereas a condensate 
would be a string in a background of p particles. The bound state improves stability 
for a much larger range of parameters than the condensate. Their conclusion was that 
the bound state definitely improves the stability of the Z string, but they could not 
find a stable solution within physical bounds (the stable solutions had too high a value 
for Weinberg angle). But they suggest a more extensive exploration of parameter space 
and an increase of the charge on the string. Since the fields involved in our context 
are to the order of fermion masses, it may be that these numbers improve things even 
more. For dynamical H and G, we have a situation corresponding to this bound state, 
and where the leptoquarks are non-dynamical, we have the condensate. 

In stability calculations, the starting point is to calculation the second differential of 
the potential with respect to the field. This is the same as the coefficient in the mass 
term and it is a significant term in dE the change in energy stored in the defect core as 
the string solution is perturbed. We saw earlier that |pp gave a positive contribution 
to the energy change by fiattening off the hump, that is, if p = 0, as in the standard 
model, the extremum at the symmetric phase is a maximum and any perturbations 
about that are obviously unstable. Whereas the presence of p changes that. It is clear 
that p ^ raises the energy in the string core since it tries to make the Higgs mass 
go positive. We conclude that the presence of p significantly improves the stability of 
the Electroweak string. We can of course repeat the same argument for case [D) and 
indeed (-D2), in which the physical parameters involve the degrees of freedom in G, 
which are related to the fermion masses. To determine exactly how significant these 
findings are, and to write a quantitative statement about the stability, which depends 
crucially on the values of the physical parameters involved, some numerical analysis 
would be required including solving a Schrodinger equation. 
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5.4 The "Higgs force" and torsion gravity 



Connes creates a picture of two identical space-time sheets, one labelled left and the 
other right. The geodesic distance between the two sheets is given by the metric (the 
Dirac operator) on the noncommutative internal space. (See for example [FCG] and 
|GM] ). When this metric is flat, the sheets are a constant distance apart. There is 
a connection one-form describing parallel transport between the two sheets. As the 
metric is fiat, it is a flat connection, its curvature is zero. The 'flat' Dirac operator 
of internal space is then 'fluctuated' by the noncommutative algebra, (this is called 
an 'internal fluctuation') and by the equivalence principle, the condition that the 
curvature is zero is relaxed, and it takes on the role of the Higgs field. In this way, the 
Higgs scalar degrees of freedom are picked up. The Higgs is promoted to a dynamical 
variable with non-zero curvature. So, now the metric is no longer flat and the geodesic 
distance between the two space-time sheets is determined by the eigenvalues of the 
fluctuated Dirac operator, which is of course, the fermion mass matrix. In the standard 
model, the internal space does not consist only of two points left and right, but models 
the whole of charge space. The mass matrix has a non-degenerate set of eigenvalues, 
as for each fermion there is a different distance between left and right. 

The finite spectral action, (the Higgs potential HI), in which when varied with respect 
to all degrees of freedom in Dp, one gets a fully degenerate set of fermion masses, may 
be written as: 

Sv(mggs) = Tt{MM* - If 

whereas if the equations of motion had a solution such that all the fermion masses 
were different, it would take a form of the type: 

Scurvature = Tt{MM* + GG* - if. 

In other words, there would be additional mixing terms between a set of new fields 
and the Higgs. Inspecting the derivation of the action in the spectral action principle, 
we may interpret this term as the curvature term for internal space. So, to get a set 
of different masses, we need extra curvature giving rise to extra Higgs force via the 
equivalence principle. The action [H] indeed provides this type of additional terms, but 
G is only non-zero in the vacuum of M if the leptoquarks are non-propagating, that 
is, if they do not possess internal Einstein's field equations of their own. 

Since the action is asking for an extra term, we are led to twist the Dirac operator 
one more time and use the Weintzenbock equation such that the sought for additional 
curvature term may appear but we do not wish to add another G field as in 15.41 in 
an ad hoc manner. At least the leptoquarks have a raison d'etre. We explain below 
why the extra term could just as well be due to torsion gravity instead of gravity 
from Riemannian geometry so that the new internal space gravity force is more like 
torsion gravity than curvature gravity. Arguably, it makes sense to interpret these 
extra gravity-making terms as torsion rather than curvature because: 
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1. G satisfies a definition of an internal space torsion-type parallel transport sug- 
gested below, 

2. the terms are intimately connected with the Electroweak string solutions, which 
are analogous to line defects in laboratory experiments that can source torsion, 
and 

3. to include torsion in the formalism is another natural extension, not an "add- 
on." It is the removal of the constraint that unnaturally forces internal space 
torsion to be zero. 



The commutative spectral triple is a Riemannian geometry, which means that the con- 
nection is defined to have zero antisymmetric part or torsion. The phrase Riemannian 
geometry is often taken to be synonymous with general relativity and the Higgs force 
is an analogy of gravity. A Riemann-Cartan geometry is one in which the connection 
is allowed to have a non-zero antisymmetric part and in this way it is a generalisation 
of general relativity. The gravitational field strength as usual, is described by the 
connection and it has not only a curvature 2-form but also a torsion 2-form. Torsion 
happens when we parallel transport two infinitesimal vectors in the tangent space at 
a given point along one another and the resulting parallelogram does not close. The 
torsion tensor is a measure of the non-closure of the resulting parallelogram. (This 
contrasts the notion of curvature, which characterises the difference between a vector 
before and after it is parallel transported around a loop on a curved manifold.) In 
"torsion gravity" theories, the source of gravity is torsion instead of curvature. Unlike 
in torsion gravity theories, instead of intrinsic spin as the source of torsion, we have, 
if it be stable, the Electroweak string. So, by including torsion in the 0-dimensional 
spectral triple of the model we are studying, we are not making an unnatural addi- 
tion, but simply removing the ad hoc constraint that the torsion be zero. The article 
|WRj shows that the space-time torsion in noncommutative geometry does not have 
the status of dynamical field as its kinetic term disappears (wherein there is a product 
with the symmetric metric tensor and the antisymmetric torsion). 

M and M* are parallel transports between left and right and MM* is a parallel 
transport transcribing a complete loop, from left to right and back to left. It appears 
in the internal space curvature term identified in the action, if is a map from right 
to right, and after a Wick rotation (as partially described in[2]), G and G* are maps 
between left and left or right and right. So H and G are parallel transports in the 
tangent space and the action terms in G*G should be a measure of torsion. 



5.5 Vacuum manifold metric and string energetics 

In |EVGj it was found that the string solutions define circular paths (geodesies) on 
the vacuum manifold^ We have illustrated this point earlier in this section. For this 

^^In |EVG| sphalerons (Electroweak monopoles) are also studied. 
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reason, Kibble and Lepora expected to find that the spectrum and properties of the 
strings should correspond to the geometry of the vacuum manifold. (The metric on 
the vacuum manifold is given by the Euclidean inner product at each tangent space.) 
So, the boundary conditions, which specify the solutions themselves, correspond to 
the geodesies with respect to that metric. They found two inequivalent metrics, one 
corresponding to each of the gauge and scalar sectors. 

What we find is that since there is a different metric on the internal space for each 
fermion (see the first paragraph of I5.4p . which in turn determines the metric on the 
vacuum manifold, there must be a set of geodesic paths with respect to each of those 
metrics. (The three-sphere + |02p = 1 becomes |ap + |6p = 1 — \g\'^ — \h\'^ for the 
up quark and so on. To avoid repetition of these equations, we refer to |TSMj for the 
corresponding metric pertaining to the other fermions.) For a given metrical structure, 
one has a corresponding curvature of the vacuum manifold to which the energy of the 
embedded vortex is associated. Given a stable string, these metrics and geodesies 
paths correspond to that string's boundary conditions. In other words, since there 
is a discrete set of masses, there is a discrete set of metrics on the vacuum manifold 
with respect to which geodesic paths on the vacuum manifold correspond to embedded 
string solutions in this extension to the Weinberg-Salam theory, whose stabilities are 
connected to the fermion masses as we have seen. 
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